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Abstract 

We present the studies of the decays B — > ai(1260) 7r and ai(1260) K within the framework 
of QCD factorization. Due to the G-parity, unlike the vector meson, the chiral-odd two- 
parton light-cone distribution amplitudes of the a\ are antisymmetric under the exchange of 
quark and anti-quark momentum fractions in the SU(2) limit. The branching ratios for a\ tt 
modes are sensitive to tree-penguin interference. The resultant B(B° — ► afn T ) are in good 
agreement with the data. However, using the current Cabibbo-Kobayashi-Maskawa angles, 
(3 = 22.0° and 7 = 59.0°, our results for the mixing-induced parameter S and a e ^ differ 
from the measurements of the time-dependent CP asymmetries in the decay B° — > afn T 
at about the 3.7a level. This puzzle may be resolved by using a larger 7 > 80°. For 
a\K modes, the annihilation topologies give sizable contributions and are sensitive to the 
first Gegenbauer moment of the leading-twist tensor (chiral-odd) distribution amplitude of 
the a\ meson. The B — > a\K amplitudes resemble the corresponding B — > ttK ones very 
much. Taking the ratios of corresponding CP- averaged a\K and ttK branching ratios, we can 
extract information relevant to the electroweak penguins and annihilations. The existence 
of new-physics in the electroweak penguin sector and final state interactions during decays 
can thus be explored. 
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I. INTRODUCTION 



The first charmless hadronic B decay involving a l 3 Pi axial-vector meson that has been 
observed is B° -> af (1260) tt t 0, 0, B H8 , 

which goes through b — > uud. The measurements 
of time-dependent CP asymmetries in hadronic B decays originating from b — > uud can 
provide the information directly related to the Cabibbo-Kobayashi-Maskawa (CKM) weak 
phase a = arg(— V t dV t * b /V uc iV* b ) (or called 2 ), for which some results have been given from 
the data of B — > 7r + 7r~, p ± n T and p ± p =F 1. The BaBar collaboration recently reported the 
observation of B° — > af (1260)7r T , including CP violating parameters, branching fractions, 
and a e ff, where the bound on the difference Act = a — can be constrained by using the 
broken SU(3) flavor symmetry 0, S] or isospin analysis ja, 10, 11 . 

In this paper, we present the phenomenological studies of B — > a\Ti and a\K within the 
framework of QCD factorization, where the former processes are tree-dominated, while the 
latter are penguin-dominated. The ai(1260), which will be denoted by a x for simplicity, 
is the l 3 Pi state. Due to the G-parity, the chiral-even two-parton light-cone distribution 
amplitudes (LCDAs) of the a\ are symmetric under the exchange of quark and anti-quark 
momentum fractions in the SU(2) limit, whereas, unlike the vector meson, the chiral-odd 
two-parton LCDAs are antisymmetric. Ref. 12] is the only literature so far for the calculation 
of LCDAs of l 3 Pi axial-vector mesons. The large first Gegenbauer moment of the leading- 
twist tensor distribution amplitude of the a\ meson [l2| could have a sizable impact on the 
annihilation amplitudes. On the other hand, it is interesting to note that, for the axial- vector 
mesons with quantum number l 1 Pi, their chiral-even LCDAs are anti-symmetric under the 



exchange of quark and anti-quark momentum 



odd two-parton LCDAs are symmetric [12j, Il3|. The hadronic B decays involving such a 



tactions in the SU(3) limit, while the chiral- 



meson are sensitive to the new-physics search |14 

Because the axial-vector and pseudoscalar penguin contributions interfere constructively 
in the dominant decay amplitudes of B — > aiK, for which the K is emitted and a\ shares the 
same spectator quark within the B meson, the B —>■ a±K amplitudes resemble very much the 
corresponding B — > nK ones. Moreover, larger CP asymmetries could be expected in the 
a\K~ and a\K modes due to the much lager decay constant of the ax(1260), as compared 
with 7r.fr channels. 

To resolve the puzzle about the observations of the decays B — > ttK and within the 
Standard Model (SM) joljjjsoHie approaches were proposed, including considerations of final 
state interactions (FSIs) |16|, ll?J, Ll8|] , and use of SU(3) flavor symmetry to extract hadronic 
parameters from the titt data and then to predict Kn channels 19|, |20|, |21[ . On the other 



hand, it was argued that new-physics with a large CP-violating phase may exist in the 
electroweak penguin sector 19, 



20 



22 



The present studies for B — > ain and a\K modes 
can offer further tests for the above theories. 

The layout of the present paper is as follows. In Sec. II, we discuss light-cone distribution 
amplitudes for an axial-vector meson. A brief description for applying QCD factorization to 
the decays B — > a\K and a\K is given in Sec. Ill, where some relevant formulas are collected 



in Appendices |A] and [Bl In terms of the notations and (3f , which were given in Ref. 23 
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one can find that the amplitudes for AP modes have the same expressions with those for 
PP and VP modes (where A = the axial-vector meson, P = the pseudoscalar meson, and 
A = the vector meson). Sec. IV contains the numerical analysis for the branching ratios 
and CP asymmetries. Our conclusions are summarized in Sec. V. 



II. TWO-PARTON LCDAS OF THE a ± AND PROJECTION OPERATORS ON 
THE LIGHT-CONE 

For decays involving an axial-vector meson (denoted as A) in the final state, the QCD 
corrections can turn the local quark- ant iquark operators into a series of nonlocal operators 
as 



(A(P,X)\q la (y)q 2S (x)\0) 

» 



j fi ( ( f*W ? 

- / due^ p y +Upx nf A m A i ^75 H u )+ <l X) 159^\u) 

A Jo [ \ pz 

/ ) 8a 

where the chiral-even LCDAs are given by 

(A(P, XMiyh^MlO) = ifAm A \ due^™^ \p, + g ( f\u) , 

j o I pz J 

(2.2) 

(A(P, XMivhM^M = -ifAm A e, upa e* ( ^z° £ duS""**** (2.3) 

with u (u — 1 — u) being the momentum fraction carried by qi(q~2), and the chiral-odd LCDAs 
are given by 

(A(P, A)|9!(y)(V7552(aO|0> = /l [ due**"**** {(e^p, - c£> P/t ) $ ± («), 

+ m \^y Z iP^-P^) hf{u)^ (2.4) 

(p)/ 



(A(P,A)|g 1 (y)75g 2 (x)|0) = /im^( A )^ / due ^vy+n px ) ^±1 , ( 2 . 5 ) 

j o I 



Here, throughout the present discussion, we define z = y — x with z 2 = 0, and introduce the 
light-like vector p^ = P^ — m^z^/ (2Pz) with the meson's momentum P 2 = m 2 A . Moreover, 
the meson polarization vector e* has been decomposed into longitudinal (e| M ) and transverse 
projections defined as 



= ~Pz 1 ^ ~ ~Pz Z ^ ) ' e J-M — S _ c Im ' ( 2 -6) 
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respectively. The LCDAs $||, $j_ are of twist-2, and g±\g±\ h± , hup of twist-3. For the a\ 



meson, due to G-parity, ^\\-,g^ and gf 1 are symmetric with the replacement of u <-> 1 
whereas hfp and hup are antisymmetric in the SU(2) limit 12]. Here, we restrict ourselves 
to two-parton LCDAs with twist-3 accuracy. 

Assuming that the axial-vector meson moves along the negative z-axis, the derivation 
for the light-cone projection operator of an axial-vector meson in the momentum space is 
in complete analogy to the case of the vector meson. We separate the longitudinal and 
transverse parts for the projection operator: 



» 



u. 



M. 



See 



5a- 



(2.7) 



where only the longitudinal part is relevant in the present study and given by 



Mi 



J A m A (e*n + ) 



4 

+ iE 



^_75$||(w) 



4 



2E 



- e>>75 n - n + h 



(*), 



u 



dv($±(v) -hP(v)) a^n 







u 



dk 



75" 



k=up 



with the momentum of the quark q\ in the A meson being 



k 2 



(2.8) 



(2.9) 



for which E is the energy of the axial-vector meson and the term proportional to k\ is 
negligible. Here, for simplicity, we introduce two light-like vectors nl = (1,0,0,-1), and 
n+ = (1,0,0,1). In general, the QCD factorization amplitudes can be written in terms of 
the form J* du Tr(M^ - • •) . 

In the following, we will give a brief summary for LCDAs of the a\ mesons, for which 
the detailed properties can be found in Ref. [12!]. ^\(u) can be expanded in Gegenbauer 
polynomials: 



$ll ±) (u) = 6uu 



i=0 



|| (±),a lrl 3/2 



1 C- /z (2n- 1) 



(2.10) 



For the $ 



ai 

IK-L) 



due to the G-parity, only terms with even (odd) Gegenbauer moments 



survive in the SU(2) limit. In the present work, we consider the approximations: 



u 



6uu \ 1 + al' ai - 



u 



18a^' ai uu(2u 



5(2m-1) 2 -1 
1). 



(2.11) 
(2.12) 



Note that we have defined = f ai since the product /^ai" ,ai always appears together. 
Neglecting the three-parton distributions and terms proportional to the light quark masses, 
we can relate the twist-3 distribution amplitudes to the twist-2 ones by Wandzura-Wilczek 
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relations 12, 24 and then obtain: 



hf(v) = (2v-l) 



"^du-f^du 



hp(v) 



u 



u 



<li! — I — - du 



u 



du(§_i_(u) - h\i(u)) = vv 



v u 



it 



The normalization conditions for LCDAs are 



du$\\(u) = 1, 



= (2v-l)$ (v) 



C?M — 



-2$ a (^; 

v U 



du 



= VV$Jv) 



f 1 du®±(u) = 0. 
Jo 



(2.13) 



(2.14) 



duhf\u) = 0, 



duh^(u) = 



(2.15) 



For the pseudoscalar meson (P) with the four-momentum P M , the light-cone projection 
operator in the momentum space reads 



M p = i^E ^„7 5 $p(m) + 



0» 



• eA ix d 4> p {u) 



(2.16) 



k=up 



where «p = m 2 P /(mi + m 2 ) is proportional to the chiral condensate (with mi 2 the masses 
of quarks) and the approximate forms of LCDAs that we use are 



$ P (u) = 6mm 1 + 3of (2m - 1) + a{ 



5(2w- l) 2 - 1 



$ p (u) 



w(l — w). 



(2.17) 



III. DECAY AMPLITUDES 



Within the framework of QCD factorization, in general the effective weak Hamiltonian 
matrix elements for B — > M\M 2 decays can be expressed in the form 23 



i—i G f 



(M 1 M 2 \H eS \B) = -^J2 X p(MiM 2 \T a p+T b p\B) 



(3.1) 



where X p = V p bV* q with q = d or s, M 2 is the emitted meson, and M 1 shares the same 
spectator quark within the B meson. Considering a generic 6-quark decay, T^ p describe con- 
tributions from naive factorization, vertex corrections, penguin contractions and spectator 
scattering, whereas Tg p contain the weak annihilation topologies. 
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For B decay processes, the QCD factorization approach advocated in |25|, |26] allows us to 
compute the nonfactorizable corrections in the heavy quark limit since only hard interactions 
between the (BM\) system and M 2 survive in the mj — > 00 limit. Naive factorization is 
recovered in the heavy quark limit and to the zeroth order of QCD corrections. In this 
approach, the LCDAs play an essential role. In the present study using the notations af 



and f3f given in Ref. 23j, the amplitudes for AP modes have the same expressions with 
those for PP and VP modes; B — > a\7T, a\K decay amplitudes in terms of af and (3f can 
be obtained from B — > pre, pK [23J by setting p — > a\. However, one should note that the 
determination of the relative signs of the detailed amplitudes behind the coefficients af and 
(3f is non-trivial. 



A. Decay amplitudes due to T A % 



In general, Tj^ can be expressed in terms of ca^(M 1 M 2 ) X^ BMl ' M2 \ where c contains 
factors of ±1 and ±l/\/2 arising from flavor structures of final-state mesons, aij are functions 
of the Wilson coefficients (see Eq. ( 13. 7h ). and 



X (SA, P) = { p {q) \ {V -A),\0){A(p)\(V-Ar\B{p B )) 

= -2if P m A V BA (q 2 )(e* {x) p B ) , 

x (bp,A) = ( A (q)\(V - A),\0)(P(p)\(V - Ay\B(p B )) 

= -2if A m A F BP (q 2 )(e\ x) p B ). 



(3.2) 



(3.3) 



Here the decay constants of the pseudoscalar meson P and the axial-vector meson A are 
defined by [27J 



CP(p)|g27„75<Zi|0> = -ifppp, (A(p, A)|g 27M 7 5 gi|0) = if A e£>. 



(3.4) 



The form factors for the B — > A and P transitions are defined as 27 

2 



(A( Pl X)\A,\B(p B )) = 1 



m B + m A 



^u a ^ x) p B fA BA (q 2 



(A(p,\)\VJB(p B )) 



(m B + m^Vf^q 2 ) - (e^>p B )(p B + p),-^- 



+2m A - 



(P(p)\V,\B(p B )) 



Q 

(Pb + p)u ~ 



m r — m. 



F 1 BP (q 2 ) + 



Tn B + m A 



\F BP {q% 



where q = p B - p, V BA (0) = V BA (0), F BP (0) = F BP (0) and 

yBA[q2) = ™B + m A vBA{q2) _ ™*-™V *V)- 



2m A 



2m A 



(3.5) 



(3.6) 
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The coefficients of the flavor operators a P can be expressed in terms of af as follows: 



ai(MiM 2 ) 
a 2 (M 1 M 2 ) 
a&{M x M 2 ) 

c^(MiM 2 ) 



3,EW 



4.EW 



(MxM 2 
(M 1 M 2 



a l (M 1 M 2 ), 
a 2 (MxM 2 ) , 

(MxM 2 ) - al{M 1 M 2 ) , 
j a2(AfiAf 2 )+r^og(MiAf 2 ) 

a^MxM,) - a p 7 {M x M 2 ) , 

f a? (M 1 M 2 )+rf^af(M 1 M 2 
1 a? (M 1 M 2 )-rfaf(M 1 M 2 



for MiM 2 = 
for MxM 2 = 

for M X M 2 
for MiM 2 



AP, 
PA, 

= AP. 
= PA. 



(3.7) 



where 



2m 2 P 



x m 6 (/i) 



m b (jj)(m 2 + mi)(ji) 
2rriA 



(3.1 



The effective parameters af in Eq. (13.71) to next-to-leading order in a s can be expressed in 
forms of 23 . 



a P (MiM 2 ) 



Ci + 



Ci±l 



c i± i C F a 



Ni{M 2 



An 



Vi(M 2 ) + —H^M.M^ + Pf(M 2 ), 



(3.9) 



where q are the Wilson coefficients, p = u,c, Cf = (N% — l)/(2N c ) with N c = 3, the upper 
(lower) signs refer to odd (even) i, M 2 is the emitted meson, Mi shares the same spectator 
quark within the B meson, and 



for i = 6, 8, and M 2 

1 for the rest. 



ai, 



(3.10) 



Vi(M 2 ) account for vertex corrections, Hi(M x M 2 ) for hard spectator interactions with a hard 
gluon exchange between the emitted meson and the spectator quark of the B meson and 
Pi(M 2 ) for penguin contractions. The detailed results for the above quantities are collected 
in Appendix |A] Note that in the present case, some relative signs change in Hi as compared 
with the PP and VP modes. 



B. Decay amplitudes due to 7g p — annihilation topologies 



The B — > AP amplitudes governed by the annihilation topologies read 

G G r ^ 

—7= \(AP\T B P \B) = -%—=f B f A f P Ap & i b i + e 5^3,EW + e 6 6 4 ,Ew 

V * p=u,c V ^ p=u,c L i=l 



(3.11) 
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where the coefficients are process-dependent and weak annihilation contributions are pa- 
rameterized as 



b 2 

&3.EW 
bdEW 



N? 



ciA\, 



Cf 
N 2 

a 



h = j^ [c,A\ + c 6 4] 



c 9 A\ + c 7 (Ai + A f 3 ) + N c c 8 Ai 



c 10 A\ + c s Al\ . 



(3.12) 



The subscripts 1,2 and 3 of A% denote the annihilation amplitudes induced from (V— A)(V— 
A), {V — A)(V + A) and (S — P)(S + P) operators, respectively, and the superscripts i and 
/ refer to gluon emission from the initial and final-state quarks, respectively. For decays 
B — > AP, the detailed expressions for A 1 ^ are given in Appendix [B] j3f(M 1 M 2 ) are defined 
as 

-ifBfu 1 fM 2 



/?f(M 1 M 2 



X(BM lt M2) 



IV. NUMERICAL RESULTS 



A. Input parameters 

In the numerical analysis, we use the next-to-leading Wilson coefficients in the naive 
dimensional regularization (NDR) scheme 281 ] . The relevant parameters are summarized in 
Table fl] |29|, |3Cj, [3l|, [32j, [331] . The value of fs that we use is consistent with the lattice average 
341 ] . The current value of F Bn (0) becomes a little smaller, and is more suitable to explain 
the 7T7T data jfjj]. We use the light-cone sum rule results for the B — » it, K 31] and B — > a\ 



32| transition form factors, for which the momentum dependence is parametrized as 35 



/0? 2 ) = /(o) 



i 



_J- r BZ(Y) q 2 /m 2 B * 

q 2 jm\, 1 - a B z(Y)q 2 /m B 



(4.1) 



where m B * is the lowest-resonance in the corresponding channel. Note that since the mass 
of the a\ meson is not small, we have, for instance, [F 1 B7r (m^ i )/i 71 ] S7r (0)] 2 ~ 1.2. It means 
that the q 2 dependence of B — > tt, K form factors cannot be ignored in the prediction. As 
for the B — > a± form factor, its q 2 dependence can be negligible due to the small mass of 
pseudoscalar mesons. However, to be consistency, I also consider its q 2 dependence in the 
analysis. Our light-cone sum rule result for V Bai (0) is a little larger than the previous QCD 
sum rule calculation, 0.23 ± 0.05 36(. It is interesting to compare with other quark model 



calculations in the literature. The mag nitude of V£ a H0) is about 0.13 and 1.02 ~ 1.22 in the 



quark model calculations in Ref . [37j] and Refs. [38|, ]39] , respectively. The magnitude of the 
former is too small and the latter is too large if using them to compute the branching ratios 
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of B — > afn T and then comparing with the data. The values of the Gegenbauer moments 



of leading-twist LCDAs for the a\ meson are quoted from Ref. 121 ] . The integral of the B 



meson wave function is parameterized as 25 



1 - p A B 

where 1 — p is the momentum fraction carried by the light spectator quark in the B meson. 
Here we use A B (1 GeV) = (350 ± 100) MeV. 

There are three independent renormalization scales for describing the decay amplitudes. 
The corresponding scale will be specified as follows: (i) the scale p v = m b /2 for loop diagrams 
contributing to the vertex and penguin contributions to the hard-scattering kernels, (ii) pu = 
\J p v Ah for hard spectator scattering, and (iii) pa = V p v Ah for the annihilation with the 
hadronic scale Ah ~ 500 MeV. We follow j25| to parameterize the endpoint divergences Xa = 
Jq 1 dx/x and X H = Jq 1 dx/x in the annihilation and hard-spectator diagrams, respectively, as 

X A{H) = In (^) (1 + PA(H)e^), (4.3) 

with the unknown real parameters Pa, Ph and 4>a,4>h- We adopt the moderate value 
Pa,h < 0.5 and arbitrary strong phases 4>a,h with pam = by default, i.e., we assign a 
50% uncertainty to the default value of Xa(h) (with pa,h = 0) jio], 41]; with the allowed 



ranges of pa,h, the theoretical predictions for nK modes are consistent with the data. Note 
that the a±K rates could be sensitive to the magnitude of pA- 

B. Results 

We follow the standard convention for the direct CP asymmetry 

- B(g' -. S) - B[B° -, /) 

The branching ratios given in the present paper are CP-averaged and simply denoted by 
B(B — ► /). The numerical results for CP-averaged branching ratios and direct CP asymme- 
tries are summarized in Tables [Til and IHH respectively. The results for time-dependent CP 
parameters of the decay B(t) — > afir^ are shown in Table HVl 



1. B — > a\-K 

The decay of the B° meson to afir^ was recently measured by the BaBar and Belle groups 
A recent updated result by BaBar yields [4] 

B(B° -> afn T -> 7T T 7r ± 7r ± 7r T ) = (16.6 ± 1.9 ± 1.5) x 10" 6 . (4.5) 
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Assuming that B(af — > n T 7r ± 7r ± ) equals to B(af — > 7r T 7r 7r°) and B(af — > (Stt^) equals to 
100%, they have obtained 



B{B° 



(33.2 ±3.8 ±3.0) x 10" 



(4.6) 



Very recently, the measurements of time-dependent CP asymmetries in the decay B — > afn 



have been reported by the BaBar collaboration [5j. From the measurements, the individual 
branching ratios of B — > afir~ and a^7r + can be obtained. As given in Table [Til our 
theoretical results are in good agreement with experiment. It was shown in Ref. |42j that 
three-parton Fock states of M 2 can give non-small corrections to o^, so that \a^\ — 0.30, If so, 
we can expect B(B — > a^ir ) ^ 1.6 x 10~ 6 , which can be tested in the future measurement. 

The B — > a\TT amplitudes are analogous to the corresponding B — > pn ones [43|. The 
tree(T)-penguin(P) interference depends on the sign of sin 7 (where V u b = \ V u b\e~ n ) and the 
relative sign between Re (a?) and Re^f); for sin 7 > 0, the T-P interference is destructive 



in B — > afn'* 1 , B 



a?7r" 



", while it is constructive in B 



Because the amplitudes 



of a\iT and pir modes are dominated by the terms with a\ and af, and ~Re[a p A (Ttai / 



Re[a4(ai7r)/3] 
relations, 



Refa^vrp)] « — Re[«4(p7r)] ~ —0.034, one can easily obtain the following 



B(W 



B(B 
B(B° 



B(B° 
B[B- 



B[B- 
B(B~ 



B(B- 



Qi 7T 



p 71 



p^ix 



P°7T 



P~7T 



< 



< 



> 



Ff*(m*)f p . 



[mi 



Vr 



Bai , 



, F^(ml 1 )f ai \'B(B°^aT7r+) 



F^(ml)f p 



B(B l 



which can offer constraints on the magnitudes of f ai and VJ- 



Bai 



m 



B(B° -> a^7x+)/B{R 



B(B" -> aj7T+) 
B(B° -> at ir-) 



a{ it 

^B-n/^2 



IS 



(4.7) 



p 71 ' j 

Moreover, the ratio 



+ 2 



' V Bai (ml)U 



1 ±Re 



1 Re 



X t ( a 4 (7rai) - a 4 (ai7r) ± (3 3 (7rai) - /3 3 (ai7r) 
An \ ai(7rai) 
A(vrai) 



a\(irai 



+ 0(a 



v 

4,EW) 



$4> /^3,EW) ^4,ew) > (4-8) 



which is not only sensitive to the form factor and decay constant of the a\ meson but also 
to the weak phase 7. The measurement of the above ratio allow us to obtain the further 
constraint on the value of 7. 

The large direct CP asymmetries may result from the non-zero value of the weak annihila- 
tion parameter (p^) and its corresponding phase. See Table HTT1 With default parameters, the 
direct CP asymmetries for at tt~, a^7r + , a^Ti , a\-K~ are only at a few percent level, whereas 
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TABLE I: Summary of input parameters. 



Running quark masses [GeV] and the strong coupling constant [23, 291 



m c (m c ) 



m s (2GeV) 



{m u + m d )/{2m s 



a s (l GeV) 



1.3 



0.09 ±0.01 



0.0413 



0.497 



Wolfenstein parameters for the CKM matrix elements [30] 



A 



X 



P 



V 



0.806 



0.2272 



0.195 



0.326 



Decay constants for mesons [MeV] [12, 29 . 401 



f. 



K 



f. 



B 



fai 



131 



160 



195 ± 10 



238 ± 10 



Form factors and parameters for their q dependence [31, 321 



Ff n (0) = 0.26 ±0.03 



ubz = 0.40 



r BZ = 0.64 



m\ = tub* = 5.32 GeV 



Ff K (0) = 0.33 ± 0.04 oiBZ = 0.95 r B z = 0.52 m x = m B * = 5.41 GeV 



Vf ai (0) = 0.28 ± 0.03 



a Y = 0.90 



r Y = 0.65 



m i = tub* = 5.32 GeV 



Gegenbauer moments for leading-twist LCDAs of mesons at scale 1 GeV [12, 33\ 



a. 



K 



,K 



1*1 



a: 



Ml 



a 



Ml 



0.25 ±0.15 



0.06 ± 0.03 



1.05 ±0.15 



-0.03 ±0.02 



-1.04 ±0.34 



it can be very remarkable for the a°7r° mode. At the present time, the large errors in the 
measurements for CP asymmetries do not allow us to draw any particular conclusion in 
comparison with theoretical predictions. (See Tables HTll and HVl) 



2. Time- dependent CP for B(t) — > a l ir 



±<rrT 



Following Ref. [8|, we define 

A + = A(B° -> at*-) , A_ = A(B° -> a^7i + ) , 

A + = A(B° -> a^vr+) , A_ = A(B° -> a+Tr") . (4.9) 

Neglecting CP violation in the B° — B° mixing and the width difference in the two B° 
mass eigenstates, time- dependent decay rates for initially B° decaying into af tt t can be 
parameterized by 



T(B°(t) - af^) = e - Tt ± (\A ± \ 2 + \A T \ 



x 1 + (C ± AC) cos Amt - (S ± AS) sin Amt\ , (4.10) 

where 

I 4 J2 _ \~A |2 

C±AC= ! % ; _ T ' , (4.11) 
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TABLE II: CP-averaged branching fractions for the decays B — ► ai(1260)7r and oi(1260)i^ (in units 
of 1CT 6 ). The theoretical errors correspond to the uncertainties due to variation of (i) Gegenbauer 
moments, decay constants, (ii) quark masses, form factors, and (iii) \b,Pa,h, 4>a,Hi respectively, 
added in quadrature. 

Mode Theory Expt. (BaBar) [4,5] Expt. (Belle) [3] 

S° - a+vr- 8.7t 2 2-ii-ii 12.2 ±4.5 

S°-_;> + 25.1±l;t±l;l±?:l 21.0 ±5.4 

S° -> a^vrF 33.8±|-|±|-|±f-^ 33.2 ± 5.0 48.6 ± 5.6 

R- _> „0_- 7 O+0.3+1. 7+1.3 

rs — > a 1 7T ' -°-0.3-1.5-0.9 

R° _> „0-^0 fi n +0.4+5.6+9.7 

D > aji\ °- u -0.4-2.6-3.1 



d— . _— i^" i Q -1 +1.3+15.5+24.5 

-D ^ «i A 19 -1-1.3- 7.8-11.0 

R- _> r fiK~ 11 a+1-0+8.7+13.1 

£> > Oxii ll-8_i.o-4.6- 4.8 



TABLE III: Direct CP asymmetries for the decays B -> ^(1260)^ and a x (1260) K (in %). See 



Table iLtl for errors. 



Mode 


Theory 


BaBar [5, 6] 


Mode 




Theory 


If -> a+vr- 
5° af^ 
B° -> a?vr° 
£?- -> a^vr 
J3- -> a?vr- 


n+0.1+0.3+20.1 
°- z -0.0-0.5-19.5 

1 7 +0.1+0.1+13.6 

-0.1-0.0-13.4 
RQ q+5.4+6.9+25.0 
Dy -°-6.1-8. 9-74.7 
n 4 +0.4+0.2+ll.l 
u -^-0. 4-0.1-11.1 
n r+0. 5+1. 5+13.0 
u -°-0.3-2.4-14.6 


7 ±21 ±15 
15 ± 15 ± 7 


iTU 
£f° — ► 
£T -> 
B~ -> 


a$__- 


9 7+0.2+0.9+11.8 
z -' -0.2-0.8-11.9 
7 q+0.7+2. 1+7.6 
'• y -0. 7-2. 2-8.3 

n 7+0.0+0.1+0.6 

u -' -0.0-0.1-0.1 
o q+0.5+1. 5+12.1 
°-°-0.5-1.7-13.4 



and _ 

S ± AS. ^^ff . (4,2) 

Here Am denotes the neutral B mass difference and V is the average B° width. For an initial 
B° the signs of the cos Ami and sinAmt terms are reversed. The four decay modes define 
five asymmetries: C, S, AC, AS, and the overall CP violating A^p, 



= /LT -rL ; rL , V, • (4.13) 



1^4+ 


2 ± 


A_ 


2 _ 


A. 


2 _ 


^4+ 


2 


1^4+ 


2 ± 


A_ 


2 ± 


A_ 


2 ± 


^4+ 


2 



Two a-related phases can be defined by 

o4 = i arg(e- 2 ^ ± ^) , (4.14) 
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TABLE IV: Parameters of the time-dependent B — * afir^ decay rate asymmetries. 
S and AS are computed for j3 = 22.0°, corresponding to sin(2/3) = 0.695, and 
7 = 59.0°. See Table ILT1 for errors. 



Theory Experiment (BaBar) [5] 



o-oi Too ±S:go±8:gi -°- 07 ± °- 07 ± °- 02 

C 0.02 ±g-gg ±g'88 ±g'il -0.10 ± 0.15 ± 0.09 

S -0.551^+0-04+o.os 0.37 ±0.21 ±0.07 

AC 0.48±8-8t±o'o2±ooi °- 26 ± °- 15 ± °- 07 

-o.oi±g;88±8:88±8:8i -o.i4±o.2i±o.oe 

hntc i +0.3+0.9+4.4n,o 
_o.3-0.5-2.4J 
Q +0.6+3. 2+6. 4\o 
l lio ' y -0.6-2.1-3.6/ 

a eS (109.5±g:l±;i±l:g)° (78.6 ±7.3)° 



which coincide with a in the limit of vanishing penguin amplitudes. The average of and 
a~ s is called a e fj: 



_ a eff + a cS 1 

Q!pff = = — 

2 4 



./ S + AS \ . ( S-AS 

arcsm . + arcsm 



1-(C + AC)V \Jl-(C-ACf 



(4.15) 



The numerical results for the time-dependent CP parameters are collected in Table IIV1 The 
magnitudes of Afyp, C and AS are small in the QCD factorization calculation, where C is 
sensitive to the annihilations and can be ~ 10% in magnitude. AC describes the asymmetry 
between B(B° -> a+vr-) + B(B° -> a{n + ) and B(B° -> ar7r+) + i3(fi° -»■ a+TT"), and 
thus can be read directly from Tables [Til and IIHI Neglecting penguin contributions, S and 
a C H, which depend on a(— n — (3 — 7), coincides with sin 2a and a, respectively, in the 
SM. Using a = 99.0°, i.e., 7 = 59.0°, the numerical results for S and a e g differ from the 
experimental values at about the 3.7a level. This puzzle may be resolved by using a smaller 
a = Ti — (3 — 7 < 78°. In Fig. [T], we plot S versus 7 (and a), where we parameterize 
V ub = 0.00368 e'^. The best fitted value is 7 = (87±f )°, corresponding to a = (71±l 3 )°, for 
P = 22°. 



3. B -> _i(1260)__ decoys 

The decays _? — > oi-fT are penguin-dominated. Because the dominant axial-vector and 
pseudoscalar penguin coefficients, 04(01 K) and ag(aiK), are constructive in the a\K modes, 
B — > a\K and the corresponding B — > ttK decays should have similar rates. It is instructive 
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a 

158° 108° 58° 8° 
1.5 




0° 50° 100° 150° 
7 



a 



158° 108 c 



58 c 



8° 



150° 



% 100° 




50° 100° 150 c 



r 



FIG. 1: S and a e s versus 7 (and a) for adopting [5 = 22°. The solid curves are obtained by using 
the central values (default values) of input parameters. The region between two dashed lines is the 
theoretical variation within the allowed range of input parameters. 



to consider the four ratios: 



Ri 



R-2 



Ra 



B{B -> aJK- 



B(B° 



Ti+K- 



Bai , 



m 



K) 



a 



i( ai K)\- 



\ a%{ ai K) otfaK) ) 

B(B~ -> a^Tf) 
B{B~ 



7T 



K 



m 



KJ 



a,%(aiK) 



k""kj / \ al(nK) J 

1 + 2 R J Pj( ai K) + P^jaJC) _ Pj(nK) + P^ttK) ' 

\ ot%{a\K) a%(wK) 



B(B° -> agiT 
£(5° - 



al(aiK) \ 
a c A (nK) ) 



1 -3 Re 



a. 



3,EW 



(tfai) 



a 



a%{a\K) 

1 tfC 



-n 



3,EW 



(Kit) 



"''2 



B(B- 
B[B- 



7T 



0#- 
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vo Bal (m 2 K 

F n B7r (m 



+2 Re 




a1{a\K) 



a%(nK) 



1 + 3 Re 



3,EW 



3,EW 



g(7Tg) + P^AttK) 

a c 4 (nK) 




-r 2 



(4.16) 



where 



ri = 



r 2 = 



ai 



/ai 



Vo Bai (ml)f K 



F BK (m 



Ft {ml) f K 



i 1.9, 
1.1, 



(4.17) 
(4.18) 



and the dots stand for the neglected terms which are numerically estimated to be less than 
1% in magnitude. The ratios #1,2,3,4, which are very insensitive to 7, are approximately 
proportional to [V Bai (m 2 K ) / (F Bv (m 2 K )} 2 and receive corrections mainly from the electroweak 
penguin and annihilation topologies. The value of the annihilation /3 3 is sensitive to ai ,ai . 
The contributions originating from electroweak penguin and annihilation amplitudes can be 
further explored by taking into account the following measurements for ratios, 





Ri 




R~2 ~ 


Ri 


R3 


R 2 




Ri 


Ra 


R 3 


~ R~2 : 



1 -3 Re 



Pi EW ( ai K) Pl EW (nK) 



a%(aiK) 



_L £ ^ 6 R e 



— i ^ — t ^ 1 + 3 Re 



«3,Ew(^ fl l) «3,EW 
■ = T\ ■ 

a%(aiK) 
a 



a%(itK) 
(Kir) 



,Ew( Ka i) 



a%(nK) 



-r 2 



(4.19) 
(4.20) 



2\R 2 #4 




3,Ew(^) r 



(4.21) 



Although the above ratios are parameterized according to the QCD factorization, they can 
be treated in a model-independent way. It is worth stressing that because ®±(u) is anti- 
symmetric under interchange of the quark and antiquark momentum fractions in the SU(2) 
limit, the weak annihilations (and hard spectator interactions), which could contribute siz- 
able corrections to the decay amplitudes, enter the B — > a\K amplitude in a very different 
pattern compared with B — > ttK decays. More relevant information about and a l ,ai can 
thus be provided by the measurement of R\j ' R 2 . 

With default parameters, the direct CP asymmetries are analogous to the correspond- 
ing B — > irK modes; because A C ps are dominated by Re(V t * d V tb ) lm.(a% + ^3) Im(V^* d (V^,) 
times Re[«i + a 2 F BK f a J (V Bai f K )} and -Re[a 2 F BK f a J(V Bai f K )] terms for a^K~ and 
a\K modes, respectively, their direct CP asymmetries are thus a little larger than the 
corresponding ttK modes in magnitude due to the decay constant enhancement. Note that 
the value of (3% is sensitive to the first Gegenbauer moment of ^°^(u) and the annihilation 
parameters pa and 0^. On the other hand, an outstanding problem is the determination 
of the signs for direct CP observations in the irK modes. The experimental results are 
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A C p(tf -> 7T+K-) = -0.095 ± 0.013 and A CP (B- -»■ vr ^-) = 0.046 ± 0.026 |6j. Some 
proposals, for instance the contribution due to new-physics in the SM electroweak penguin 
sector 0, 0, Q or due to FSIs HI 181 1 . were advocated for the resolution. The ratio mea- 
surements for R1/R2 — R3/R4, R1/R3, and -R4/-R2 directly probe the electroweak penguins. 
Moreover, the approximate relation given in Eq. (14.211) will be violated if the FSI patterns 
are different between a\K and irK modes. 



V. CONCLUSIONS 

We have studied B — ► ai(1260) 7T, ai(1260)X decays. This paper is the first one in the 
literature using the QCD factorization approach to study B — > AP decays. Interestingly, 
due to the G-parity, the leading- twist LCD A &± of the ai(1260) defined by the nonlocal 
tensor current is antisymmetric under the exchange of quark and anti-quark momentum 
fractions in the SU(2) limit, whereas the ^jj 1 defined by the nonlocal axial-vector current is 
symmetric. The large magnitude of the first Gegenbauer moment (a^ ,ai ) of $^ could have 
a sizable impact on the annihilation amplitudes. If one ignores i.e., letting a^' ai = 0, 
with default parameters (where pA = 0), the branching ratio for a®K° mode becomes 1.8 
times smaller, while the changes of branching ratios for aiir and the remaining a\K modes 
are at the level of 5% and 10%, respectively. 

Our main results are summarized as follows. 

• Our results for B(B° — > a^vr - , a^7r + ) are in good agreement with the data. Theoret- 
ically, the rates for B — > ai(1260) it are close to the corresponding ones for B — > pw. 
The differences between the above two modes are mainly caused by different mag- 
nitudes of form factors (V Bai and Aq P ) and decay constants (f ai and f p ), and by 
different patterns of tree-penguin interference. For sin 7 > 0, the T-P interference is 
destructive in B —>■ a : f7T ± ,B~ — > a\n~ , but constructive in B~ — > a^ir . Because 
the amplitudes of aiir and pre modes are dominated by terms with ot\ and ctf, and 
Re[c^(7rai)] ps Re[a^(ai7r)/3] w Re[c^(7rp)] ps -Re[c^(p7r)] ps -0.034, we obtain the 
relations as given in Eqs. fl4.T[) and (14. 8p . Thus estimates for form factors and decay 
constants as well as the weak phase 7 can thus be made from these ratio measurements. 

• For CP asymmetries, the large experimental errors do not allow us to draw any par- 
ticular conclusion in comparison with theoretical predictions. The time-dependent CP 
asymmetry measurement in B° — > afn T can lead to the accurate determination of the 
CKM angle 7. Using the current fitted value 7 = 59.0°, i.e., a = 99.0° corresponding 
to (3 — 22.0° in the SM, our results show that S and a e s differ from the present data at 
about the 3.7a level. This puzzle may be resolved by using a larger 7 > 80°. Further 
measurements can clarify this discrepancy. 

• The branching ratios for the decays B — > a\ir and a\K are highly sensitive to the 
magnitude of V^ Boi (0). Using the LC sum rule result, V Bai (0) = 0.28 ± 0.03 Q, the 
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resultant branching ratios for af ir T modes consist with the data very well. Neverthe- 
less, the value of V^ ai (0) is about 0.13 and 1.02 ~ 1.22 in the quark model calculations 
in Ref. 37 1 and Refs. (ss|, 3^], respectively. If the quark model result is used in the 



calculation, B(B° — > af ir T ) will be too small or large as compared with the data. 

The B — > a\K amplitudes resemble the corresponding B — > nK amplitudes very much. 
Taking the ratios of corresponding CP-averaged a\K and ttK branching ratios, we can 
extract information about the transition form factors, decay constants, electroweak 
penguin (ag EW (Kai)), and annihilation topology (/^ EW (aiK)). See Eqs. fl4.19p - fl4.2ip . 
Thus, the possibilities for existing new-physics in the electroweak penguin sector and 
for final state interactions during decays can be explored. 



Note added. Recently Belle has updated the following measurement 44]: B(B —>■ afn + 



a 1 7r + ) = (29.8 ± 3.2 ± 4.6) x 10 6 which is in good agreement with our result. On the 
other hand, BaBar has reported new measurements on a^Ti^ , a^ir and afK~, a~{K modes 
0, G3> where B(B° af K~) = (16.3 ± 2.9 ± 2.3) x 10" 6 is also in good agreement with 
our prediction, whereas the central values of branching ratios for the remaining modes are 
about 2 ~ 3 times larger than our predictions. The latter discrepancies should be clarified 
by the improved measurements in the future. 
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APPENDIX A: THE COEFFICIENTS a 



p 



In the below discussion, we set = $ p . In Eq. (13. 9p . the expressions for effective 
parameters a\ are 



a\{M x M 2 ) 



Ci + 



Cj±l 



Ni{M 2 



Cj±i C F a s 



N r 4tt 



47T 



Vi(M 2 ) + _^(M 1 M 2 )j + Pf (M 2 ). 



(Al) 



Ni is given in Eq. (I3.10p . The vertex corrections have the same expressions as those for VP 
modes 23|] with LCDAs of the vector meson being replaced by the corresponding ones of the 
a\ meson. For the penguin contractions Pf (M 2 ), one can perform the same replacements but 
needs to add an overall minus sign to P<f (ai) and P$(ai). Hi(MiM 2 ) have the expressions: 



H i (M 1 M 2 



Ah 



o 1 — p 
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for % = 1 - 4, 9, 10, 
H i (M 1 M 2 ) 



x / dv du 



uv 



± r 



Mi 



uv 



if B fM 1 fM 2 f 1 d ?l(P) 



X{BMx,M 2 ) J r 1 - p 



dv I du 



uv 



uv 



(A2) 



(A3) 



for i = 5,7, and H i (MiM 2 ) = for z = 6,8, where the upper (lower) signs apply when 
M 1 = P (Mi = A). Here $f (p) is one of the two LCDAs of the B meson [25 . 



APPENDIX B: THE ANNIHILATION AMPLITUDES Ah f 



For (see Eq. (13.121) ). some signs change in comparison with the results of B — > PP 
and PV . We obtain 



4 = Tia s J^ dxdy ^ 2 {x)^ 1 {y) 
A{ = 4 = 0, 



1 1 

+ 



A\ = na s I dxdy \ <$>f\x)^(y) 



y(l - xy) x 2 y 



1 1 

+ 



rfrf $ m2 (x) S mi (y) — 



x(l — xy) xy 2 



rfrf 4> m2 (x)$ mi (y) 



Al = ±na s fdxdy {rf $f (x) $ mi (y) - % + rf $f (y) $ 
jo [ ^Z/(,t xyj 



2 

xy 
2x 



2(l + x) 



4 = ±* a .£dxdy | r f $f (x)$ mi (y)^^-rf $f (y)$ m2 (x 



'"- , ~ y xy(l - xy) 

2(1 + y) 



x^y 



xy 2 



(Bl) 



where the upper (lower) signs apply when Mi = P (Mi = A) and the detailed definitions of 
the distribution amplitudes of the axial-mesons have been collected in Sec. HI] Again, here 
we have set $j| = 4> p . 

Using the asymptotic distribution amplitudes of ^jj 1 (w) and $p(w), and the approximation 
&±(u) = 18uu(2u — l)cii ,ai , we obtain the annihilation amplitudes 



A\ « 67ro; s 



A* 2 « 67ra s 



3 X.-4+ 



7T 



ai- a V^r^(X A - 3) 



3 ( A A - 4 + y ) - ai' ai r a x ^X A (X A - 3) 



.4' 



±67TQ;o 



7T 



<( Ai - 2X A + — ) + 3a^' ai r^ ( A^ - 2A A - 6 + 



7T 



4 « 6vra s (2A A -l) A A - 3a±' ai r? (X A - 3) 



(B2) 

(B3) 

(B4) 
(B5) 



18 



B. Aubert et al. [BABAR Collaboration] , |arXiv:hep-ex/0408021 



B. Aubert et al. [BaBar Collaboration], |arXiv:hep-ex/0507029| 



K. Abe et al. [Belle Collaboration], arXiv:hep-ex/0507096 



B. Aubert et al. [BABAR Collaboration], Phys. Rev. Lett. 97, 051802 (2006) 
[arXiv:hep-ex/0603050| . 



B. Aubert et al. [BABAR Collaboration], arXiv:hep-ex/0612050 



See Heavy Flavor Averaging Group, http:/ /www. slac.stanford.edu/xorg/hfag/, and references 
therein. 

M. Gronau and J. Zupan, Phys. Rev. D 70, 074031 (2004) [arXiv:hep-ph/0407002| . 

M. Gronau and J. Zupan, Phys. Rev. D 73, 057502 (2006) [arXiv:hep-ph/0512148| . 

Y. Grossman and H. R. Quinn, Phys. Rev. D 58, 017504 (1998) [arXiv:hep-ph/9712306| . 

J. Charles, Phys. Rev. D 59, 054007 (1999) |arXiv:hep-ph/9806468| . 

H. J. Lipkin, Y. Mr, H. R. Quinn and A. Snyder, Phys. Rev. D 44, 1454 (1991). 

K. C. Yang, Nucl. Phys. B 776, 187 (2007) [ arXiv:0705.0692l [hep-ph]]. 

K. C. Yang, JHEP 0510, 108 (2005) | arXiv:hep-ph/0509337] . 

K. C. Yang, Phys. Rev. D 72, 034009 (2005) [Erratum-ibid. D 72, 059901 (2005)] 
|arXiv:hep-ph/0506040| . 

P. K. Das and K. C. Yang, Phys. Rev. D 71, 094002 (2005) | ]arXiv:hep-ph/0412313| . 

W. S. Hou and K. C. Yang, Phys. Rev. Lett. 84, 4806 (2000) [Erratum-ibid. 90, 039901 (2003)] 

[arXiv:hep-ph/9911528] . 

C. K. Chua, W. S. Hou and K. C. Yang, Mod. Phys. Lett. A 18, 1763 (2003) 
|arXiv:hep-ph/0210002| . 

H. Y. Cheng, C. K. Chua and A. Soni, Phys. Rev. D 71, 014030 (2005) |arXiv:hep-ph /04093 1 7| . 
A. J. Buras, R. Fleischer, S. Recksiegel and F. Schwab, Phys. Rev. Lett. 92, 101804 (2004) 
[arXiv:hep-ph/0312259] . 

A. J. Buras, R. Fleischer, S. Recksiegel and F. Schwab, Nucl. Phys. B 697, 133 (2004) 
larXiv:hep-ph/0402112l. 



R. Fleischer, S. Recksiegel and F. Schwab, arXiv:hep-ph/0702275 
S. Baek and D. London, |arXiv:hep-ph/0701181[ 
M. Beneke and M. Neubert, Nucl. Phys. B 675, 333 (2003). 
S. Wandzura and F. Wilczek, Phys. Lett. B 72, 195 (1977). 

M. Beneke, G. Buchalla, M. Neubert, and C.T. Sachrajda, Phys. Rev. Lett. 83, 1914 (1999); 
Nucl. Phys. B 591, 313 (2000). 

M. Beneke, G. Buchalla, M. Neubert, and C.T. Sachrajda, Nucl. Phys. B 606, 245 (2001). 
M. Bauer, B. Stech and M. Wirbel, Z. Phys. C 34, 103 (1987). 

G. Buchalla, A.J. Buras, and M.E. Lautenbacher, Rev. Mod. Phys. 68, 1125 (1996). 
Particle Data Group, W. M. Yao et al, J. Phys. G 33, 1 (2006). 



CKMfitter Group, J. Charles et al, Eur. Phys. J. C 41, 1 (2005); http://ckmfitter.in2p3.fr 



19 



[31] P. Ball and R. Zwicky, Phys. Rev. D 71, 014015 (2005) [arXiv:hep-ph/0406232| . 
[32] K.C. Yang, in preparation. 

[33] P. Ball, V. M. Braun and A. Lenz, JHEP 0605, 004 (2006) [arXiv:hep-plT70 603063 1 ■ 

[34] M. Bona et al. [UTfit Collaboration], JHEP 0610, 081 (2006) |arXiv:hep-ph/0606167| . 

[35] P. Ball, Phys. Lett. B 644, 38 (2007) [arXiv:hep-ph/0611108| . 

[36] T. M. Aliev and M. Savci, Phys. Lett. B 456, 256 (1999) [arXiv:hep-ph/9901395] . 

[37] H. Y. Cheng, C. K. Chua and C. W. Hwang, Phys. Rev. D 69, 074025 (2004) 

|arXiv:hep-ph/03103"59] . 
[38] D. Scora and N. Isgur, Phys. Rev. D 52, 2783 (1995) |arXiv:hep-ph/9503486| . 
[39] A. Deandrea, R. Gatto, G. Nardulli and A. D. Polosa, Phys. Rev. D 59, 074012 (1999) 

[arXiv:hep-ph/98112"59] . 

[40] H. Y. Cheng, C. K. Chua and K. C. Yang, Phys. Rev. D 73, 014017 (2006) 

|arXiv:hep-ph/0508l"04] . 
[41] H. Y. Cheng, C. K. Chua and K. C. Yang, larXiv:0705.3079l [hep-ph]. 
[42] K. C. Yang, Phys. Rev. D 69, 054025 (2004) | ;arXiv:hep-ph/0308005] . 
[43] W. S. Hou and K. C. Yang, Phys. Rev. D 61, 073014 (2000) [arXiv:hep-p h/9908202|. 
[44] K. Abe et al. [Belle Collaboration], [arXi\^070JL3J279] [hep-ex]. 
[45] D. Walker [BABAR Collaboration], [ar2Oyiy708 i 0050J [hep-ex]. 

[46] D. Brown, invited talk presented at the XXIII International Symposium on Lepton and Photon 
Interactions at High Energy, August 13-18, 2007, Daegu, Korea. 



20 



